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Abstract 



The 4>Tr + Tr~ production in pp annihilation at rest is strongly enhanced 
by a two step mechanism with intermediate KKirir states. The relative 
yield of the <fi production due to the resonant final state interaction de- 
f"^. ' creases with increasing total energy of the pp system. 
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%*• 1 Introduction 

The production of mesons in low energy pp annihilation (see |], §, [3|, f|, || § 
and references therein) is strongly enhanced in some channels for which the 
production is expected to be suppressed on the tree level according to the Okubo- 
Zweig-Iizuka (OZI) rule |7|, ||. These conspicuously OZI rule violating reactions 
are of considerable theoretical interest since they can be used for studying reaction 
mechanisms and the nucleon structure, in particular the problem of an intrinsic 



ss component ITOL 11, |12|. Some well-known OZI-rule breaking mechanisms are 



two-step processes with ordinary hadrons in the intermediate state |13| , therefore 
their role in the pp annihilation must be investigated. 

The two step mechanisms in nucleon- ant inucleon annihilation have been stud- 



ied for various final states containing mesons [14], [15], [IB, [17], WE, FS, [20, 211 E2 



23]] . It has been found that in the 07r and 00 channels, which show the most 
dramatic violation of the OZI rule, two-step mechanisms play an important role 
and result in a cross section comparable with the experimental data. However, 
for a number of reactions, including pp — > 0p, 07T7T, 0cu, no significant contribu- 
tions from two-step mechanisms have been found. The goal of this paper is to 
reanalyze the role of two-step mechanisms for the pp annihilation into the 07T7t 
channel. To this aim we go beyond the approximations used in previous studies. 



1 



The resulting cross sections turn out to have the right order of magnitude. The 
plan of the paper is as follows. Section ^ gives a brief summary of the experimen- 
tal data and the previous calculations. A simple model illustrating the two-step 
mechanism with a resonant final state interaction is considered in Sec. |[ The 
reaction pp — ► 07T7T is discussed in Sec. |] where the two-step mechanism with 
the KKnn intermediate state is studied in detail. The summary of the results is 
given in Sec. ||. For earlier reviews concerning the two-step mechanisms we refer 

to [pi m eh, mm. 



2 The OZI Rule Violation in pp 



Table |l] gives a summary of the experimental data on the OZI rule violation in the 
07r + 7r~ and (pp channels at different energies. The tree level expectations based on 
the deviation of the uo — (p mixing angle from the ideal one are (0 — 0j) 2 ~ 4 • 10~ 3 
[34}] . The OZI rule violation in these channels is therefore rather moderate^. The 
experimental branching ratios for the reactions pp — > nir(p, p<p, ttttuj, pu, ttttKK 
at rest are given in Tab. 0. 



Table 1: The relative probabilities for the (pp, oop and (pii + ii~, ujii + ix~ channels 
in pp annihilation vs. antiproton momentum. Phase space corrected ratios of 
yields are also tabulated. 



BR(<f>TT + TT~ 



Pp(GeV/c 



measured 



BR(4>P) in3 
BR(ujp) ' 



PS corrected 



BR(u 



measured PS corrected 



1 



■ 10 3 



Ref. 



BE! 



(gas) 

(gas/LX) 

(liq.) 

(liq.) 

(3 atm) 

0.76 

1.2 

2.3 

3.6 



6.3 ±1.6 
7.5 ±2.4 



9±5 



22 ± 13 



13 ±4 



25 ± 15 



7.0 ± 1.4 
4.9 ±0.8 
5.9 ±0.9 
10.0 ± 2.4 
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+3 



'-4 

21 ±5 

h4 



9 



15 ±3 
10.3 ± 1.6 
12.5 ±2.0 

19 ±5 
19" 



h5 

30 ±7 

h5 



-7 



12 



1H 



-9 



H 

33j 

H 



Until now the importance of the two-step mechanisms in the reactions pp — > 
07T7T, pp — > (pp, pp — > cf)uj remained unclear. The main attention was focused 
on two-particle intermediate states (two meson doorway approximation). The 

1 This magnitude of the OZI rule breaking is often termed nondramatic, contrary to the 
cases where the <f> production exceeds the estimate from <fi — uj mixing by more than one order 
of magnitude. 
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Table 2: The experimental branching ratios for pp annihilation at rest into 
ttttKK, p<p, and 7171(f). 



Reaction 


BR 


Condition 


Ref. 


PP — > 7T + 7T~(f) 


4.6(9) ■ 10" 4 
5.4(10) ■ 1(T 4 
7.7(17) • 10 
4.7(11) • 10- 4 
6.6(15) • 10- 4 

3.5(4) • 10 4 

7/ c \ 1 n— 4 
o. ^OJ • 1U 


liq. 

gas 
gas La 

S 

P 
liq. 
gas 


N 

rTTTi 

ED 
E 
i 

i 

HI 


+ 1 

pp — ► 7T 7T <p^K s K L 




liq. 




VD — > /5(f) 
ff 


3.4(8) • 10" 4 


Eras 


im 




4.4(12) • 10~ 4 


gas./LX 
'So 
3 Pj 


1 




3.4(10) • 10- 4 
3.7(9) • 10~ 4 


i 

E 


— > 7T + 7T~U) 


6.6(6) • 10- 2 


liq. 


pqi 


pp — > pc<j 


5.4(6) • 10~ 2 
3.0(7) • 10- 2 
6.4(11) • 10- 2 


gas. 
S 
P 


|29| 
|29| 
1291 




2.3(2) ■ 10~ 2 


liq. 


13J 


PP — ► 7t + 7l~ KgKl 


2.41(36) • 10" 3 


liq. 


P6| 




2.26(45) ■ lO" 3 


liq. 


[35] 



contribution of the K*K+KK* intermediate state was considered in the unitarity 
approximation in |L7 1 ; the calculated branching ratios for the (pp and (puj channels 
were found to be about two orders of magnitude smaller than observed. The 
possibility of a large contribution of the pu intermediate state for the 07T7T channel 
was pointed out in [ T6 ] , however the estimate was done in unitarity approximation 
and neglecting spin. Since the unitarity approximation is likely to be suppressed 
by threshold factors, the off-mass-shell contributions can be large, but they are 
also known to be model dependent. Intermediate states with more than two 
particles should also be taken into account. 

Of special interest are the intermediate states containing KK, because the 
KK production is not OZI suppressed and the final state interaction (FSI) effects 
are strong if the KK system is produced in the region of the resonance. The 



two-step mechanism pp — > ttttKK — > 7171(f) was estimated as well in |16[ using 
the unitariry approximation and neglecting the spin structure of the amplitude, 
the resulting branching ratio being about one order of magnitude smaller than 
the experimental value. In this paper we present a more detailed calculation of 



3 



this mechanism which includes spin effects and off-shell contributions. It shows 
that KK rescattering in the KKtttt system leads to a significant enhancement of 
the 07T7T production in pp annihilation at low energies which has the right order 
of magnitude. Note that this straight KK — cf) rescattering mechanism cannot 
contribute to two body final states in the two meson doorway approximation. 



3 Resonant Rescattering Mechanism 

To illustrate the basic features of the resonant rescattering mechanism we consider 
a three-particle decay a —>■ 123 where particles 1 and 3 interact via resonance b. 
The total decay amplitude corresponds to the sum of the two diagrams shown 
in Figs.|I]a,b. As shown in Fig.^|, all rescattering terms 1 + 2 «-> b are taken into 
account by including the mass operator U b (pl) in the resonance propagator. Here 
and below p n denotes the four-momentum of particle n, pi = (pi + P2) 2 = $i 2 . 
The imaginary part of the mass operator U b (si 2 ) is determined by the width of 
the decay a — > 1 + 2: 

Imn 6 (s 12 ) = -y/si2~T b (s 12 ) , (1) 



r b (s 12 ) = (2) 

where g b is the coupling constant for the decay b —>■ 12 (for the sake of illustration 
all particles are assumed to be spinless in this section), and ^12(^12) is the three 
momentum of particles 1 and 2 in their CMS: 

^12(^12) = ^\J(s 12 - (mi + m 2 ) 2 )(l - (mi -m 2 ) 2 /s 12 ) . (3) 

If the resonance b is narrow, then the following approximation for the resonance 
propagator can be used: 

(4) 



pl - ml - U b (Pb) Pb- m l + im b^b 

where nib and r b = T b (ml) are the physical mass and the physical width of the 
resonance. 

The amplitude corresponding to the resonant rescattering diagram in Fig.|l|b 
has the form 

Qa9bPl2{Sl 2 ) 9b 



71-123 " "(< + P{S 12 )Y g Sl2 - m l + imbVb 



S12 ~ml + im b T b 



+ f3{s 12 ))— v - 2 | - (5) 



where g a is the coupling constant for the decay a — > 123, and (3(si 2 ) is the ratio 
of the real to imaginary part of the loop in the diagram Fig.lIJb. The imaginary 
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Figure 1: (a) The amplitude T°_ >123 for the decay a — > 123 in plane wave approx- 
imation, (b) The decay amplitude T^ff 123 with resonant final state interaction, 
(c) The unitarity approximation T^ 123 for the resonant FSI amplitude. 




Figure 2: The equation for the propagator of resonance b coupled to the two- 
particle channel (1 + 2). 
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part of the loop, which corresponds to the particles 1 and 2 on the mass shell 
(Fig.|I]c), is determined by the coupling constant and the phase space for the 
decay b — > 1 + 2. The real part of the loop is divergent and must be regularized, 
e.g. by subtraction or by introducing form factors in the vertices: 

9a -> 9aF a (s 12 ) , (6) 
g b -> g b F b (s 12 ) . (7) 

The loop in the diagram of Figjl|b can be calculated using the dispersion relation 

yfs^ IT J(m 1+m2 )2 S - S12 + «e 



Instead of specifying the subtraction point or the form factors one can consider 
(3 as a model parameter. Neglecting the real part {j3 = 0) corresponds to the 
unitarity approximation for T^ 123 for the resonant FSI amplitude (Figfjc). 
The total amplitude a — > 123 has the form 

T a ^ 123 = g a + T r a e J 123 = g a A(s 12 ) (9) 

where A(s 12 ) is the enhancement factor resulting from the resonant final state 
interaction: 



, (z + P)m b T b \ gig - mg - 2m b A 

S12 - mj + zm b r b / S12 - m£ + zm^rf, 



According to Eqs.(|9| JT0|) the total amplitude has a pole at s± 2 ~ (m& — ^r;,/2) 2 
and a zero at S12 ~ (m;, + A) 2 where 

A = (3T b /2 (11) 

The zero of the total amplitude results from the interference between the resonant 
term and the nonresonant background, and it is essential for providing the correct 



asymptotic behaviour of the enhancement factor at large s: A(s) s ^5° \ ( see flgg 
and references therein). 

The differential decay rate has the form 



dr o _ 12 3 = (2n) i ^-\A((p 1 +p 2 ) 2 )\ 2 d^ 3 ( Pa ,p 1 ,p 2 ,p 3 ) (12) 

where d& 3 (p a ,pi,p 2 ,ps) is the differential 3-body phase space, see Appendix A. 
The distribution in the invariant mass of the pair (1+2) is given by 

= (27r) 7 -^|A( Sl2 )| 2 $ 2 (m a , v ^ r :,m3)$ 2 ( v ^,m 1 ,m 2 ) (13) 
OjS\ 2 ^fn a 
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Here and below $ n (m,mi,...,m n ) denotes the total n-body phase space, see 
Appendix A. 

The resonance approximation for the production of the particles 1 and 2 
corresponds to the case when only the resonant term T^ff 123 in Eq.(|9|) is taken 
into account: 



rribTb ds 



^r a ->63->i23 — r a ^ b3 - j~r~- — wr? — ~ (i^) 

\ s 12 - mf, + vm b l b \ z 7r 

r n -,M = / dr. 



a^b3 ~ J a ^ b3 ~ 

q 2 

= (27r) 7 -^(7rm 6 r 6 )(l + /3 2 )$ 2 (m a , m b , m 3 )§ 2 (m b , mi, m 2 )(15) 

where (3 = /3(m 2 ). The ratio of the resonant production rate r a _,.&3 to the total 
rate of the decay a — > 123 in the plane wave approximation (without FSI) 

r°_i23 = (27r) 4 ^-$3(m a ,m 1 ,m 2 ,m 3 ) (16) 

is given by the formula 

Ta^b3 _ ^ mb Y b )(i | / ? 2| W 3 ^(m a ,mj,,m 3 )$ 2 (m 6 ,m 1 ,m 2 ) _ 



r°_i23 $ 3 (m a ,mi,m2,m 3 
= (irm b T b )(l + [3 



o, $ 2 (m a ,m 6 ,m 3 )$2(w6,wi,m 2 ) 



$ 2 (m a , m 12 , m 3 )$ 2 (mi 2 , mi, m 2 )dm\ 2 

Note that the off-shell contribution leads to the enhancement factor (1 + j3 2 ). 
The resonance approximation ( |i~4"|) can be used only for \/3\ ^> 1, otherwise one 
cannot neglect the interference of the resonant term with the background. For 
example, if /? = 0, then the peak structure arising from the resonance pole is 
completely suppressed by the zero in the amplitude (|9l,|T0D at s i2 = m\, so that 
the production cross section features a dip instead of a peak in the resonance 
region. 

The generalization to the four-particle decay a — > 1234 with final state in- 
teraction in the system (1 + 2) is straightforward. In particular, the ratio of the 
resonant production rate r a _^ 34 to the total rate of the decay a — > 1234 in the 
plane wave approximation is given by the formula 

r a ^634 / n W1 . n 2\ (27r) 3 $ 3 (m a ,m b ,m3,m 4 )$ 2 (m 6 ,mi,m 2 ) _ 



(7rm b T b )(l+(3 2 ) 



1^1234 $4(m a ,mi,m 2 ,m 3 ,m 4 
= (7rm b T b )(l + 



$ 3 (m a , m 6 , m 3 , m 4 )$ 2 (m b , mi, m 2 ) 



$ 3 (m a , mi2, m 3 , m 4 )$ 2 (mi 2 , m x , m 2 )dm\ 2 



7 



The differential decay rate has the form 



rfr a ^ 1234 = (2n) 4 7 ^\A((p 1 +p 2 ) 2 )\ 2 d^ A (p a ,p 1 ,p 2 ,p 3 ,p i ) . (19) 
2m a 

Equations (|H],|r8|) have a very simple interpretation: the rate of resonant produc- 
tion is proportional to the fraction of the total phase space which overlaps with 
the resonant region, with the extra factor (1+/3 2 ) accounting for the off-mass- shell 
effects. 

Applying Eq. ([17],|18|) to the case of pp annihilation into KKp and KKmr we 
get the following estimates for the <p production: 



(7rm^)(l+^ 



(27r) 3 $ 2 (2m p , m , m p )$ 2 (^, m K , m K ) 



Fpp^KRp ' $ 3 (2m p} m K} m K} m p ) 

= 0.053(1 + (3 2 ) (20) 
T pp _ (j)n + 7r - _ ^^^^^ ^ ^(2n) 3 ^ 3 (2m p ,m^,m w ,m w )^ 2 (m^,m K ,m K ) 



Tpp^KKx+TT- $4(2m p , m K ,m K ,m p ) 

= 0.014(1 + (3 2 ) (21) 



Equation (|21|) at (3 — agrees with the estimate given in |T6[ . 

To demonstrate the effects of resonant FSI and to explore the difference be- 
tween the full loop calculation and the unitarity approximation we use the cou- 
pled channel model described in Appendix B. In this model the loop integrals 
are regularized by a form factor in the vertex b — > 12. Using a monopole form 
factor with cutoff parameter A gives a ratio of the real to imaginary part of the 
loop (3 = A/ P\2 ■ As an example we choose the masses and resonance parameters 
corresponding to the process pp — > KKp and calculate the <pp final state (see 
below for a proper treatment of the spin factors). The results are shown in Fig.[| 



4 The (j) meson production in the reaction pp —* 
KKtt + 7t~ at rest 

In this section we consider <fr meson production in pp annihilation at rest in the 
KKnir channel. The KK system must be in a P-wave with C-parity C K % = — 1 
to form the <fi meson. The nn system can be either in the C-odd or the C-even 
state, and, if we restrict ourselves to S'-wave annihilation, the following transitions 
are possible: 

J PC = I : PpCSJ -> {KK) L = 1>P= - 1)C =-x (nit) L=even , P=+1 , c=+1 (22) 
j™ = Q -+ . pp (i So) (KK) L=ltP= . 1)C =-i (™)l=oM,p=-i,c=-i (23) 
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0.96 1 1.04 1.08 1.12 0.96 1 1.04 1.08 1.12 

M K - K (GeV) M K - K (GeV) 



Figure 3: The distribution of the invariant mass of the KK system for the reac- 
tion pp — > KKp at rest (spinless approximation). The full rescattering calculation 
is given by the solid line: (a) /3 — 1, (b) /3 — 2, the phase space distribution by 
the dotted line. The unitarity approximation (f3 = 0) is shown in (a) by the 
dashed-dotted line, the resonance term in (b) by the dashed line. All the distri- 
butions are nomalized relatively to the total rate calculated in the plane wave 
approximation. 
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The pp( 2S+1 Lj)KKiT7!' vertices with minimum number of derivatives have the 
form: 



pp(l )— >KKtttt 



T. 



pp(Q~ + )^KK-KTT 



9o£a(h8Pil4P3PlFo 



(24) 
(25) 



where g and gi are the corresponding coupling constants, ejf is the polarization 
vector of the pp( 3 Si) state, the pi are the four- momenta of the particles in the 
final state (jp\ and P2 correspond to K and K, and ps, p^ to mr, respectively). The 
vertex form factors are denoted by F\ and Fq. In the case (p5|) , the strong tttt 
interaction due to the p meson must be taken into account, using e.g. the method 
described in Sec. [| The results for the pp annihilation into the (pp channel will 
be published elsewhere. 

In the rest of this paper we focus our attention on the production from 
the triplet S-wave pp state. The recent partial wave analysis of the 4>7nr channel 
measured by OBELIX |], [37] demonstrates that the S'-wave reaction is dominated 
by the 3 Si initial state if the irir system is produced with invariant mass below 
the p resonance, therefore our calculations can be directly compared with these 
dataQ. The coupling (0) is proportional to the relative momentum of the KK 



pair, P KK = y (Pk — Pk) 2 : /2, therefore the form factor in the annihilation vertex 
is important for damping the transition strength at large Pkk- We use the 
following parametrizations: 



KK) 



_ A 2 h 

+ pJr 
A f 

WTp 1 



(a) — monopole 

(b) — dipole 



(26) 



KK) 



where A is a cut-off parameter which will be determined later. We neglect the 
final 7T7T interaction in the S'-wave, since the energy dependence of the FSI effects 
below the /o(980) resonance is known to be rather smooth (see [5S| and references 
therein). For the sake of simplicity we also neglect the FSI effects in the Kit 
systems; this assumption, however, should be removed in further studies, because 
the K* production is rather strong in the KKim channel [^, |36| . 

The 

(27) 
(28) 



The calculation of the resonant rescattering amplitude is similar to Sec| 
result for the differential decay rate is 



dT 



KKtt+tt- 
dm KK 



(27T) 



2m 



\Ai(rn KK )\ 2 W(m m rri 



KK) 



VP 



W(m pP ,m KK ) = ^ 3 (m m m KK ,m n ,m^ 2 ('rri KK ,rn K ,'m K ) 



2 Due to low acceptance of the 4>p channel the yields of the 4>ir + ir channel measured by 
OBELIX are lower than those measured by other groups (see Table ||) . 
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where m p p = 2m p , m KR = ^/s RR = 2^jm K + P\ R is the invariant mass of the 
KK system, and the enhancement factor A\ has the form 



2 ^ (i + pjsKK^m^T^ 
s K R -ml + im 4> T 4 
h{s Kk ) = 2P KR F 1 (P KR ) . (30) 



M{s KR ) = /i(sm)-/iK)— 2 , -p , (29) 



Keeping the resonance term only gives the production rate 

rV +7r - = (27r) 7 ^(7rm b r & )(l+/5(^ i? ) 2 )/ 1 (mJ)iyK P -,m^) (31) 



which can be compared with the total decay rate to the KKtt + tt channel in the 
plane wave approximation: 

po = (7rm 6 r 6 )(l +/3 )- f v - . (32) 

KKn+n- j ft {m\ R ) W (m pP , m KR )dm\ R 

For an explicit evaluation we proceed by the following steps. First, the form 
factor (|26|) is parametrized by fitting the experimental invariant mass distri- 
bution da/dm KR in the mass range outside the resonance for the reaction 
pp — > (KK) l= itttc at rest. We use the data on the reaction pp — > K s K l tt + ty^ 
[ |36|1 which selects the final states K°K° with C kR = —1 and L kR = 1. Since the 
partial wave analysis for the final state KsKltttc is not available, we must rely on 
some assumptions about the final state decomposition. One of them is the domi- 
nance of the 5-wave annihilation, which allows two possibilities: pp( 3 Si, 1 ) — > 
(Km^c^imr)^ andpp( 1 S ,0- + ) -> (K ^°) c . = _ 1 (7r7r) c - = _ 1 . Second, we 
assume that the contribution from the singlet spin state does not bring a sig- 
nificant distortion of the shape of da/dm KR outside the <fi meson region. The 
resulting cut-off parameter is A = 0.2 GeV for the monopole form factor (p6|a) 
and A = 0.4 GeV for the dipole form factor (p6|b). 

After fixing the form factor for the pp annihilation vertex, the only parameter 
which remains to be determined is the ratio of the real to imaginary part, (3. For 
the dipole form factor fl2"6|b) in the annihilation vertex, the loop in the diagram 
Fig.|l|b is finite, and using Eq.(||) one gets 

P(*kr) = , (33) 



leading to /3 — 4.3 for A = 0.4 GeV. If the monopole form factor (pqa) is used 
in the annihilation vertex, then an additional form factor is needed for the <pKK 
vertex. Introducing a monopole form factor with the cutoff parameter A^ = A 
leads again to Eq . (|33|) , and (3 = 2.0 at A = 0.2 GeV. These estimates show that 
the real part of the loop is significant even for rather soft form factors because 
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M K - K (GeV) M K - K (GeV) 



Figure 4: The distribution of the invariant mass m K K for the reaction pp( 3 Si) — > 
KKtt + tt~: (a) f3 — 1, (b) /3 — 4. The full calculation is given by the solid line, 
the resonance term only by the dashed line, the plane wave approximation by the 
dash-dotted line, the phase space distribution by the dotted line. 



■*150 



"D 




0.9 1 1.1 1.2 1.3 1.4 1.5 

M K - K (GeV) 



Figure 5: The calculated distribution of the invariant mass m K ^ for the reaction 
pp( 3 Si) — > KKtt + tt~' {(3 = 2) in comparison with the experimental data PBfl . 
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Table 3: The calculated ratio of the (p resonance peak to the total yield for the 
reaction pp^Si) — > (KK) l=1 (ittc) l=q at rest for various values of /3, Eq.(^). 





1.0 


2.0 


3.0 


4.0 


dipole form factor 


0.046 


0.083 


0.139 


0.211 


monopole form factor 


0.059 


0.105 


0.173 


0.261 



the imaginary part is suppressed by the factor corresponding to the P-wave 
of the relative motion in the KK system. 

The results for the distribution in the invariant mass of the KK system are 
shown in Fig.f| for various values of (3. As already discussed in Sec. |3], the resonant 
FSI mechanism causes a characteristic peak-dip structrure in the invariant mass 
distribution around the (p mass. This feature can be used to distinguish the 
resonant FSI mechanism from other mechanisms of (p production, including the 
intrinsic strangeness [pi], [ll) and the two-meson doorway mechanism pp — > pu — > 
07T7T |I6(|, which give only a peak in the region of the (p meson. The calculated 



mass distribution for /3 = 2 is compared with the experimental data in Fig.[|. 
A lower limit for (3 can, in principle, be obtained from the experimental data 
using Eq.flTIl) provided the position of the zero is known. This requires the 
separation of different partial waves in the distribution of the effective mass of 
the KK system, which is not presently available. However, from the fact that the 
resonance peak in the experiment |6|, |37| is observed very close to the (p mass, i.e. 
its position is not shifted due to a nearby zero, we can conclude that (3 > 2. Thus 
both theoretical estimates and the data indicate a large value of the real part of 
the loop. The calculated ratio of the (p resonance peak to the total rate of the 
(KK) l= i(titt) l=0 channel is shown in Table |3] (the peak contribution is defined 
as the integral from the KK threshold to the zero of the production amplitude, 
see Figs.|,|). 

To calculate the absolute yield of the (pint we use the measured branching 
ratio BR(pp — > KsKl7t + tt~) = 2.4 • 10~ 3 [36 assuming a fraction of 3/4 for 



to the annihilation from the 3 Si state. Considering all charge channels in the 
intermediate state doubles the branching ratio leading to our estimate BR(pp — > 
(KK)l = i7t + tt~) = 3.6- 10 -3 . Using the results for the relative yield of the (p meson 
from Table |3| gives the absolute yield to the 07r + 7r~ channel for the pp annihilation 
at rest as BR(pp( 3 Si) — > </>7r + 7r~) > 3 • 10~ 4 for (3 > 2. This result agrees well 
with the OBELIX result for the ^tttt yield in liquid BR((j)7v + 7v-) = 3.5(4) ■ 10~ 4 § 
(this value practically excludes the contribution from the (ftp channel). The other 
measurements of BR((f)7T + 7T~) shown in Table |2| give slightly higher values due to 
the (pp contribution. The ASTERIX analysis |j] gives BR(pp(S) — > (p7i + n~) = 
4.7(11) • 10~ 4 and BR{pp{ l S Q ) -> (pp) = 3.4(10) • 10" 4 . Using these branching 
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Table 4: The calculated ratio of the resonance peak to the total yield 
B R^nn) / BR((KK) l= itttt) for the reaction pp — > (KK) l=1 tttt as a function 
of the total energy. 





2 2.5 3 4 


dipole form factor, (3 = 2 
monopole form factors, (3 = 2 


0.083 0.050 0.039 0.031 
0.105 0.058 0.041 0.027 



ratios and the fact that the singlet S-wave corresponds to about one quarter of 
the total annihilation in liquid one gets BR(pp( 3 Si) — > 07r + 7r~) = BR(pp(S) — > 
</>7r+7r-) - \BR{pp{ 1 Sq) -> <f>p) = 3.8(13) • 10~ 4 , in agreement with the OBELIX 
measurement as well as with our calculations. Note, that the old bubble chamber 
measurement BR(pp(S) — ► (j)^K s K L ^ + ^) = 1-8(3) ■ 10~ 4 gives BR(pp(S) — > 
07T+7T-) = BR(pp(S) -> K s K L n + TT-)/BR((f) -> i^-fcTz) = 5.4(9) • 10" 4 in good 
agreement with the LEAR results. 

The two-step mechanism considered is one of many possible two-step pro- 
cesses. Other four-particle intermediate states to be discussed include the K*Kinr, 



KK* 7T7T, and K*K*n7r states which arguably can lead to Lipkin cancelations [13 
However, because the branching ratio for these channels is much smaller than for 
the KKtttc channel and the dispersion integral (|8|) is saturated mainly in the low 
s region due to the relatively soft form factors, we do not expect strong cancella- 
tions between the contributions from these additional four-particle states. We do 
not include the <fi — uj mixing term in order to avoid possible double counting^. 
The contribution of the pu intermediate state, which is potentially significant 



lq| , remains to be calculated, but there is no reason to expect that there is a 
cancelation between the puj and KKmr terms (their relative phase is an unknown 
parameter) . 

Using Eq.(|27D one can estimate the dependence of the relative 07T7T yield on 
the invariant mass of the pp. For the sake of simplicity we use the phase space dis- 
tribution for the 7T7T pair and the form factors for the (KK)l=i system described 
above. The results shown in Table |] demonstrate that the fraction of the total 
phase space of the KKtitt system, which is favourable for the resonant pro- 
duction, decreases with increasing total energy. This is in qualitative agreement 
with the general trend seen in OZI rule violation in pp annihilation^. 



3 The 4> — to mixing goes via a two-step mechanism with the KK, K*K, KK*, K*K* 
intermediate states. 

4 To obtain the absolute branching ratio for the c/)tttt channel from the results shown in 
Table |J, the energy dependence of the vertex pp — > {K K) jj—ynix ; which is experimentally 
unknown, must be taken into account. Because the K and K share the same s-quark line, the 
KK pairs are predominantly produced with limited relative orbital momentum, and the energy 
dependence of the relative P-wave KK production is not expected to be strong. 
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5 Conclusion 



We have developed a general formalism describing resonant final state interaction 
in a multiparticle system. This mechanism is shown to play an important role in 
nucleon-antinucleon annihilation into OZI-rule breaking channels (07T7T, d>p). The 
resonant rescattering mechanism pp — > tt + tt~KK — > 7r + 7r~0 is studied in detail. 
Off-mass- shell contributions are found to be very significant for this process, while 
the unitarity approximation is small. The interference of the resonant term with 
the nonresonant background which is essential for providing a correct analytical 
structure of the total amplitude in the case of elastic resonant rescattering, leads 
to the characteristic peak-dip structure of the invariant KK mass distribution for 
J^x = This feature can help to distinguish the resonant FSI mechanism from 
other mechanisms of d> production when high resolution partial wave analysis 
of the KK production becomes available. The <fi peak can get an additional 
enhancement due to this interference effect. The rate of the 07T7T production at 
rest from the 3 Si state is in good agreement with the experimental data |], |6|, |3Ci|j . 
Thus no unexplained OZI rule violation is required in this case — the same 
conclusion was obtained earlier for the (bit , d>d>, dry channels as well (see [F24|, E6 



27J and references therein). The d> production due to resonant FSI decreases with 
increasing total energy as the fraction of the total phase space favourable for the 
resonance formation gets smaller. This feature agrees with the general trend of 
the OZI rule violation in the nucleon-antinucleon annihilation which is decreasing 
as a function of energy. 

It would be desirable to extend the present approach to the case when resonant 
FSI effects are taken into account simultaneously in the different subsystems 
mentioned earlier in order to achieve a unified description of the production of 
the d>, p, and K* resonances. 
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A Appendix: Phase Space 

The differential n-particle phase space d^ n (p,pi,p 2 , ■ ■ ■ ,p n ) is defined by 

n n j3 

d^ n (p, Pl ,p 2 ,...,p n ) = (2tt)^(p - II (34) 

i= i i= x 

where p is the total four-momentum of the particles with four-momenta Pi = 
(Ei, pj. The phase space reduction formula has the form 

d^ n (p,p 1 ,p 2 , ...,p n ) = d$ n -i(p,p x ,p3, . . . , p n )d<& 2 (Px, Pi, P2) {^f dpi (35) 
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The total n-particle phase space for the decay a — ► 12 . . . n is 

M W ™j = /^fap^,...,^) ,P 2 =™2 • (36) 

B Appendix: Resonant FSI in a Coupled Chan- 
nel Model 

In order to describe the resonant final state interaction in a two-particle system, 
we introduce a variant of the Weisskopf-Wigner (WW) model with two channels. 
Channel 1 is the scattering channel of interest, and channel 2 has a bound state 
| b) with bare mass (the rest of the dynamics in the second channel is ignored) . 
The only interaction in the model results from the coupling between the channels. 
The T-matrix, as a function of the invariant mass squared s, is defined by the 
Lippmann-Schwinger equation 



T u T 12 \ (0 V\ ( V\(G\{s) \ (T n T 12 
T 21 T 22 ) \V + 0J + \V + 0J\ G° 2 (s)J\T 21 T 22 



(37) 



where V is the interaction between channels 1 and 2 and G?(s) and G?(s) are the 
free Green functions: 

CfM = ^ • (39) 
s — rriQ 

Here \k) denotes the free two-particle state with relative momentum k, both 
particles have the same mass m, and s = A(k 2 + m 2 ). We assume the channel 
coupling to have the following form (all particles are spinless) 

(k\V\b) = g^k) = (40) 

where g is the coupling constant and \i characterizes the range of interaction. 
The solution for the scattering amplitude in channel 1 has the form 

/(*) = -(k\T u (s)\k) = 92 \ {kf (41) 

s — rriQ — ll(s) 

where the mass operator Tl(s) is given by 

n W = ^| G ; (S )|0 = ^JL-^ . (42) 

From Eqs. (|4l|j42]) the poles of the amplitude can be easily found. For the sake of 
illustration we consider the case of weak coupling (the bare state b is assumed to 
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be above the threshold of the scattering channel), then the scattering amplitude 
has a resonance pole at s = (nib — iTb/2) 2 m m 2 , + IT^m 2 ,). The resonance shift 
5mb = nib — ttio and width Tb due to the coupling to the open channel have the 
form 

Re n(m ) 

5m b « — (43) 

2m Q 

r ^ Im IT(m ) _ 2g 2 k 2 e(h) 



m m 

where fc& = Jm%/A~ — m 2 . 

Now we consider the decay a — > 123 described by a pointlike vertex with a 
coupling constant g a . Using standard results from scattering theory we find that 
the decay amplitude with the final state interaction between particles 1 and 2 
taken into account is proportional to the scattering wave function at zero distance 
r 12 between particles 1 and 2. The result for the above described model has the 
form 

^0-123 = ^12 = 01^) = 

= (?a(ri 2 = 0|A; 12 ) + (r 12 = 0|Go(s 12 )r( Sl2 )|A; 12 ) (45) 

\ S12 - mg - II(si 2 ) J 

where the relative momentum fc 12 is related to the the square of the invariant 
mass of the subsystem (1+2): 

s 12 = 4(m 2 + A; 2 2 ) . (47) 

Formula fl4lf) can be rewritten explicitly showing the pole-zero structure of 
the amplitude: 

s 12 - m 2 - A(s 12 ) , . 

iW 23 - 9a- — 2 ^7-^ (48) 

A <-> - Cw^y ■ (49) 



Note that the nominator in Q4"5| ) is a real function for real Si 2 which has a zero 
at S12 = m 2 « m 2 , + A(m 2 ). The distance between the zero and the pole in the 
weak coupling limit is 

m z — m b = —— . (50) 
2k h K J 
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